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ABSTRACT 

Context. Spectroscopic analysis remains the most common method to derive masses of massive stars, the most fundamental stellar 
parameter. While binary orbits and stellar pulsations can provide much sharper constraints on the stellar mass, these methods are only 
rarely applicable to massive stars. Unfortunately, spectroscopic masses of massive stars heavily depend on the detailed physics of 
model atmospheres. 

Aims. We demonstrate the impact of a consistent treatment of the radiative pressure on inferred gravities and spectroscopic masses of 
massive stars. Specifically, we investigate the contribution of line and continuum transitions to the photospheric radiative pressure. We 
further explore the effect of model parameters, e.g., abundances, on the deduced spectroscopic mass. Lastly, we compare our results 
with the plane-parallel TLUSTY code, commonly used for the analysis of massive stars with photospheric spectra. 

Methods. We calculate a small set of O-star models with the Potsdam Wolf-Rayet (PoWR) code using different approaches for the 
quasi-hydrostatic part. These models allow us to quantify the effect of accounting for the radiative pressure consistently. We further 
use PoWR models to show how the Doppler widths of line profiles and abundances of elements such as iron affect the radiative 
pressure, and. as a consequence, the derived spectroscopic masses. 

Results. Our study implies that errors on the order of a factor of two in the inferred spectroscopic mass are to be expected when 
neglecting the contribution of line and continuum transitions to the radiative acceleration in the photosphere. Usage of implausible 
microturbulent velocities, or the neglect of important opacity sources such as Fe, may result in errors of approximately 50% in the 
spectroscopic mass. A comparison with TLUSTY model atmospheres reveals a very good agreement with PoWR at the limit of low 
mass-loss rates. 

Key words. Stars: early-type - Stars: mass loss - Stars: winds, outflows - Stars: atmospheres - Stars: massive - Stars: fundamental 
parameters 


1. Introduction 


The initial mass of a star determines its evolutionary path, and is 
thus considered one of the most fundamental stellar parameters. 
Yet stellar masses derived for massive stars via spectral analy¬ 
ses are generally prone to large uncertainties, greatly hamper¬ 
ing an accurate calibration of stellar masses with their spectral 
types and evolutionary status. The so-called “mass discrepancy” 
problem, which arises when comparing stellar masses obtained 
from spectroscopy, to orbital, wind, and evolutionary models, 
has been of concern to stellar physicists for a few decades (Her- 


rero et al.|1992[|Repolust et al.|2004| Massey et al.|2012] i. While 


recent studies suggest a solution of this problem over the years 
( jWeidner & Vink||2010[ Markova & Puls|[2014) >, discrepancies 
still exist, especially in the range of giants to supergiants. 

In principle, orbital masses are independent of stellar atmo¬ 
sphere models and are therefore considered to be more robust 
(e.g., Torres et aL]|20lT) . However, orbital masses are only at¬ 
tainable in the rare case of binary systems with well-constrained 
inclination, usually owing to eclipses. Immense progress has also 
been made in the field of asteroseismology, which allows the 
measurement of stellar masses with very high accuracies from 
observed stellar pulsations. Unfortunately, the high variability in 
the outer layers of massive stars make the study of their pulsa- 
tional behavior very difficult, often hindering an effective imple¬ 


mentation of asteroseismological methods to massive stars (see 
recent review by Aerts] |20T4l and references therein). Indeed, 
spectroscopy remains the primary method to infer stellar masses 
for the majority of the massive stars. 

The spectroscopic mass of a star is derived from its radius 
R „ and surface gravity g, via M, = G 1 g, R 2 , and therefore any 
uncertainties in R t and g, propagate into uncertainties in M,. Ex¬ 
cept for the rare case where the angular diameter of a star can be 
directly measured, R* is derived from the luminosity L and effec¬ 
tive temperature T e ff of the star via the Stefan-Boltzmann relation 
R * = (47tctsb) 112 L 1 ^ 2 T- 2 . Since the effective temperature can 
in principle be constrained with decent accuracy, the main cause 
for uncertainty in the stellar radius is the error in the distance d , 
which propagates in the error in L, according to R, oc vz oc d. 
If the distance is well constrained and thus the spectroscopic ra¬ 
dius is known, the stellar mass M, is directly proportional to g. t 
and thus all uncertainties in the gravity determination propagate 
directly into mass uncertainties. 

The gravity g, is by no means a directly measurable quan¬ 
tity. It is derived by comparing synthetic spectra from model at¬ 
mospheres with observations. The surface gravity of a star de¬ 
termines the stratification of its atmospheric pressure, and most 
prominently affects the pressure-broadened wings of hydrogen 
and helium lines. The determination of g, thus relies on both the 
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pressure broadening theory adopted, as well as the model atmo¬ 
sphere. 

Different broadening theories can lead to systematic differ¬ 
ences of ~ 0.15 dex in inferred log g t values, which alone cor¬ 
responds to ~ 40% error in the inferred spectroscopic mass. It 
is an even harder task, however, to constrain the systematic er¬ 
rors that arise because of different assumptions and techniques in 
stellar atmosphere codes. The radiative pressure in massive stars 
depends on the stellar parameters as well as the opacity, which 
in turn depends on the elements, abundances, atomic data, and 
line Doppler widths. As we illustrate in this study, g * is highly 
model-dependent, and systematic errors can easily occur if the 
radiative pressure is not fully and consistently accounted for. 

Two domains can be distinguished in the atmosphere of a 
massive star: a hydrostatic domain, where gravity is balanced by 
pressure (e.g., gas pressure, radiation pressure), and a wind do¬ 
main, where the outward pressure exceeds gravity and the matter 
is accelerated. 

The spectra of O- and B-type stars are mostly formed in the 
outer layers of their quasi-hydrostatic domains. For the model¬ 
ing of such stars, a detailed treatment of the hydrostatic regime is 
imperative. Specifically, the contribution of line and continuum 
transitions to the total radiative pressure in the quasi-hydrostatic 
domain is far from negligible. Proper knowledge of the velocity 
held in the layers close to the stellar surface is a key ingredient 
for a better understanding of a variety of theoretical and observa- 


tional phenomena (see, e.g., Hamann 198l||Cidale & Ringuelet 

1993 

Owocki & Puls 1999; Cantiello et al. 2009; Shenar et al. 

2014 

Ir 

tion c 
code 

i. 

this paper, we thoroughly document the current calcula- 
f the radiative pressure in the Potsdam Wolf-Rayet (PoWR) 
(see, e.g. Grafener et akj2002; Hamann & Grafener]2003) 


and illustrate the importance of accounting for it consistently. 
While we focus on the PoWR code, the major concepts and 
equations are representative of the majority of current state-of- 
the-art stellar atmosphere codes, as we illustrate in a brief com¬ 
parison. We further discuss and quantify the large impact of a 
proper hydrostatic treatment on inferred stellar parameters, and 
particularly on the stellar mass M„. We demonstrate the sensitiv¬ 
ity of the radiative pressure to the Doppler width of spectral lines 
and to the iron abundance. Lastly, we compare our results with 
the plane-parallel TLUSTY O-grid models from |Lanz & Hubenyj 
( 2003| >, which are widely used for the analysis of hot stars with 
negligible winds. 

The structure of this paper is as follows: In Sect.[2] we briefly 
summarize the main assumptions of the PoWR code and thor¬ 
oughly discuss the treatment of the quasi-hydrostatic domain. 
The outcome of our test calculations are shown and discussed 
in Sect.[3] In Sect.[4] we compare our spectra with correspond¬ 
ing TLUSTY models, before drawing the general conclusions in 
Sect. [5] 


2. The PoWR code 

2.1. The Basics 

The Potsdam Wolf-Rayet models describe atmospheres of spher¬ 
ically symmetric stars with a stationary outflowjj To achieve a 
consistent solution, the equations of statistical equilibrium and 
radiative transfer are iteratively solved to yield the population 
numbers without the approximation of a local thermodynamic 

1 For Wolf-Rayet stars, model grids are available online at 
http://www.astro.physik.uni-potsdam.de/PoWR/ 
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equilibrium (non-LTE). The radiative transfer is solved in the co¬ 
moving frame, which avoids simplifications such as the Sobolev 
approximation. After an atmosphere model is converged, the 
synthetic spectrum is calculated via a formal integration along 
emerging rays. Some description of the PoWR code can be 
found in jGrafener et al.| ( |2002| l and |Hamann & Grafenerj ( [2004) 1 . 
The temperature stratification is updated iteratively to ensure en¬ 
ergy conservation in the expanding atmosphere, as described in 
Hamann & Grafener ( |2003j l. Recently, the PoWR code has been 
extended by the so-called thermal balance-method, which goes 
back to ideas of|Hummer & Seaton| (| 1 963]> and [Hummer (1963) 


and is described in detail for stellar atmospheres by Kubat et al. 
\\999) and [Kubatj(2001 j). This method provides better numerical 
stability in optically thin domains. 

In the comoving frame calculations during the non-LTE iter¬ 
ation, we assume that the line profiles are Gaussians with a con¬ 
stant Doppler broadening velocity which approximately ac¬ 
counts for the thermal and turbulent velocity. A constant broad¬ 
ening velocity is a well-established simplification in comoving 
frame methods. While PoWR uses directly a velocity u^op as in¬ 
put, the “Comoving Frame General” (CMFGEN) code requires 
three input parameters 7d op , Ad op , and vj - which are then com¬ 


bined to a constant velocity iq op = y2 


Martins et al. 


2002 


T'dop 

f^H^dop 


• Vj (see, e.g., 


CMFGEN description). In the FASTWINE^ 
code a velocity similar to PoWR has to be given and is referred 
to just as mircroturbulence i> tur b qPuls et al. |2005] l. 

The value of t>d op is chosen such that it approximately reflects 
the order of the averaged thermal speed combined with the mi¬ 
croturbulence, i.e.. 


^dop 




1 + y turb- 


(1) 


For O and B stars typical values for i>d op range between 10 and 
30km/s. Only for stars without photospheric lines, such as clas¬ 
sical Wolf-Rayet stars, higher values can be chosen to speed up 
the calculations without changing the emergent spectrum. The 
influence of Ud op on an O-star model spectrum is discussed and 
illustrated in Sect. 13. 41 

Pressure broadening is neglected during the comoving frame 
calculations, which is sufficient for the current applications of 
the PoWR models. For stars with considerably higher values of 
logg*, e.g. , subdwarfs, specific codes, such as the Tubingen 
Model-Atmosphere Package (TMAP) (e.g., Werner et al.|2003) 
exist which include this effect. Some codes, such as the TLUSTY 
code, have the option to switch on pressure broadening in the it¬ 
eration if needed. In the formal integration in PoWR, detailed 
thermal, microturbulent and pressure broadening are accounted 
for in a depth-dependent manner. 

The basic parameters of a PoWR model are the stellar tem¬ 
perature T *, luminosity L, mass-loss rate M, surface gravity g* 
and the chemical abundances. The stellar temperature T t is de¬ 
fined as the effective temperature of a star with the luminosity L 
and radius R, (referred to as the “stellar radius”), defined at the 
Rosseland continuum optical depth r max = 20. The total Rosse- 
land optical depth r Poss (/Lj including lines is larger. To ensure 
T max = 20 at the inner boundary, the velocity i) mm = v(R t ) is it¬ 
eratively adjusted. The surface gravity g, is defined at the stellar 
radius R, via g, = G M, r R t 2 . For OB stars, the difference be¬ 
tween the radius where r max = 20 and the “photospheric radius” 
at tr oss = 2/3 is usually very small. However, in supergiants the 
effective temperature can differ up to ~ 1 kK between these two 


acronym for “fast analysis of stellar atmospheres with winds” 
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points, and this difference in definition is apparent when com¬ 
pared with other studies. 

The density stratification in the quasi-hydrostatic domain fol¬ 
lows from an integration of the hydrostatic equation, thoroughly 
discussed in Sect. |2.2| In the wind domain, the radial wind veloc¬ 
ity v(r) is usually prescribed in the model by a so-called /3-law 


v(r) = 




( 2 ) 


where Uoc is the terminal velocity of the wind, and ft is a free 
input parameter whose value typically ranges between ft = 0.6 

With the mass-loss rate M 


and ft = 2.0 ,(e.g., Puls et al. 2008 


specified, the density stratification p(r) in the wind follows from 
the continuity equation 

M = An r 2 v(r)p(r). (3) 

In the calculation, we use complex model atoms, with a su¬ 


perlevel approach for iron group elements (see Grafener et al. 


|2002| for details), and an explicit set of quantum levels for all 
other elements. The detailed chemical composition for our cal¬ 
culations is given in Sect.[3] together with the rest of the model 
parameters. 

We do not use any clumping in the models and instead as¬ 
sume a smooth wind. The potential existence of clumping in the 
subsonic photosphere is a constant debate in the massive star 


community (see, e.g., Runacres & Owocki 2002 

Oskinova et al. 

2007; Cantiello et al. 2009; Sundqvist & Owoc 

ki|2013 1 . While 


clumping would significantly affect the derived absolute stellar 
parameters in any case, the goal of this work is to demonstrate 
effects that are independent of any clumping formalism. There¬ 
fore, we refrain from assuming a particular clumping approxi¬ 
mation and calculate our test models (see Sect. [3]! with a smooth 
wind. 


2.2. The quasi-hydrostatic domain 


The Potsdam Wolf-Rayet (PoWR) model atmosphere code was 
originally developed for WR stars where the emergent spectrum 
is formed almost exclusively in the stellar wind. For these ob¬ 
jects, it is sufficient to treat the quasi-hydrostatic domain of WR 
stars with a simple barometric formula using a constant scale 
height 


H c := 


KT, 

r 


+ V 


2 

turb 


geff R* 


(4) 


in units of R t while ensuring a smooth transition of the velocity 
field and its gradient between the two quasi-hydrostatic domain 
and the wind. Here, is the gravity corrected for radiative pres¬ 
sure in the hydrostatic domain (see details below), p the mean 
particle mass (including electrons) in units of the hydrogen atom 
mass ;;jh, and u tur b the turbulent velocity, which is a free input 
parameter. The parameter R denotes the specific gas constant for 
hydrogen, i.e., R = k^ftn^. 

For a proper treatment of O- and B-star atmospheres, the 
barometric formula with a constant scale height is not an ap¬ 
propriate solution of the hydrostatic equation. The hydrostatic 
equation is one of the fundamental equations of stellar structure. 
For massive stars, one has to account for the outward radiative 
force acting against gravitation, yielding the radial stratification 
of pressure P(r) in the hydrostatic domain. 


d P 
dr 


= -p(r)s(r)[l-r(r)]. 


(5) 


Here, p(r) is the mass density, g(r) = G M*/r 2 the gravity, and 
F( r) the ratio between the outward radiative acceleration and 
gravitational acceleration. The term [1 - I Yr)| describes the ef¬ 
fective reduction of the gravity due radiative pressure, as dis¬ 
cussed below. 

With the assumption of an ideal gas, the pressure P can be 
expressed by 


P(r) = p(r) 


RT(r) 2 

+ V. 


turb 


p(r) 

= PO) [a 2 (r) + i> t 2 rb ], 


( 6 ) 

(7) 


where T is the electron temperature. In the second line we further 
introduce the sound speed 


a{r) : = 


RT(r) 
I P(r) 


( 8 ) 


in order to simplify the expression. The turbulent velocity f,’ lurb 
is a free depth-independent input parameter reflecting a possible 
microturbulence. While in the formal integral the microturbu¬ 
lence is combined with the actual thermal velocity of each ele¬ 
ment to obtain the precise depth-dependent Doppler broadening 
velocity, it is not directly connected to the value of r dop used in 
the comoving-frame calculation. To avoid physically inconsis¬ 
tent situations, the value of a dop should be higher than for i; turh . 
So far this has been ensured by the user, but we are planning 
to implement a more detailed treatment of microtubulence in the 
comoving-frame calculations, also allowing for depth-dependent 
changes. 

The current implementation allows us to compare our OB- 
type models with those of other stellar atmosphere codes, as 
they adopt similar approaches (cf. Table [l). Stellar atmosphere 
analyses with several codes (see, e.g., |Massey et al.|2013| ) have 
demonstrated that microturbulent velocities of the order of 10 
to 20 km/s help to reproduce the observed spectral lines in cer¬ 
tain parameter regimes. It is an ongoing debate whether such 
velocities represent a real turbulent motion in the photosphere, 
as originally suggested for hot stars by Struve & Elvey ( 1934| 
and prominently reintroduced by [Hubeny et al. ( 1991] ), or if they 
are rather a “fudge factor”. Related to this open question is the 
discussion whether such a turbulent term should be included 
in the hydrostatic equation or not. In the PoWR code the term 
is included, but because of its currently depth-independent im¬ 
plementation there are no additional derivatives occuring and it 
merely leads to an offset of the sound speed. 

For simplicity, we will set fj tlllb = Okm/s in the following 
calculations, but the full result can always be recovered by re¬ 
placing a 2 with a 1 + u 2 urb . Upon dividing Eq. |5j) by -p, the term 
on the left side, which we will refer to as 


-^press(r) .— j 


p dr 


(9) 


describes the outward acceleration due to gas pressure and tur¬ 
bulent motion. Thus the hydrostatic equation can be written as 


a P ressW=£(r)[l-r(r)], (10) 

illustrating that in the hydrostatic domain the outward acceler¬ 
ation due to gas pressure has to balance gravity reduced by the 
radiative acceleration. In a strictly hydrostatic domain, we have 
therefore no net velocity, i.e., no stellar wind. 
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For an accurate description of an expanding stellar atmo¬ 
sphere, i.e., with v 4- 0, one would have to use the hydrodynamic 
equation, which can be written for a stationary, symmetric out¬ 
flow in the following form: 


1 dP dv 

+u- = -g(r)[l-r(r)]. 
p dr dr 


(id 


This means there is only one additional term in Eq. (Ill com¬ 
pared to the hydrostatic Eq. (J5J). This inertia term a m ech(r) := 
is of fundamental importance in the outer wind, but becomes 
negligible quickly below the sonic point. This also holds for a fi- 
type velocity law with f3 > 0.5 where a mcc h approaches zero for 
r —> R t . Thus in the subsonic domain the hydrodynamic equa- 
transitions into the hydrostatic form ([5]> and we obtain 


tion (11 


a quasi-hydrostatic stratification. 

In such a quasi-hydrostatic situation, we can still have a non¬ 
vanishing velocity, even though its value is subsonic and, as we 
will see in the comparison with purely hydrostatic models in 
Sect.[4j negligible for the observed spectrum. Furthermore this 
small velocity still fulfills the equation of continuity 0- This can 
be used to obtain the consistent solution for the velocity field in 
the quasi-hydrostatic domain. First, we replace the pressure gra¬ 
dient in ogress with the help of Eq. (0, 


1 d P 
p dr 

dcr a 2 dp 

dr p dr 


( 12 ) 

(13) 


We then eliminate the density in the second term by using 
the equation of continuity 0 and thus write 


da 2 

"~d7 

da 2 


2 2 

• a r v- 


1 

dr \ r 2 v 


2 a 2 a 2 dv 

- + - +- . 

dr r v dr 


(14) 


(15) 


Finally, by combining Eqs. ( fT0] > and ( |~i3j ), we obtain an equa¬ 
tion for the velocity gradient in the quasi-hydrostatic regime. 


dr 

cE 


^(l-F(r))-^ + ^ 

r- r dr 


(16) 


In principle, given a value v(R „) = y rmn at the inner boundary, 
one can obtain the velocity field in the quasi-hydrostatic part via 
direct integration of this equation. In practice, we transform this 
equation and split off the main exponential trend as this turned 
out to work better in terms of numerical stability. Thus our solu¬ 
tion for the quasi-hydrostatic velocity field is 


, a 2 (r) R; Ir-R. 

V\Jj ~ ^min o/r* \ o 


a 2 (R t ) r 2 


H r 


b{r) , 


(17) 


with a yet to be determined function b(r). We now calculate the 
derivative of Eq. (17 1 with respect to r. 


dr 2r „ r da 2 

T" = — + + u 

dr r a A dr 


1 


d b 

d7 


combine this with Eq. ( fl6| ), and obtain 

d /2 GM 
cE 

1 d /2 _ GM 
H c dr a 2 r 2 


2 1 da 2 1 

- -(- —— - + 

r a~ dr H c 


1 dcr 


2 2 (l-r(r))-- + - —. 

a A r A r a L a r 


(l-F(r)) 


(18) 

(19) 

( 20 ) 
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The righthand side of Eq. |20| i is now defined as as it 

is analogous to the definition of H c (Eq.[4]). Thus the numerical 
integration of the velocity gradient is replaced by the integration 
of 


db _ 1 1 

cb 1 ” K ~ H(P)' 


( 21 ) 


The required boundary value b(R *) = 0 follows from the inner 
boundary, where it is required in Eq. ( |T7j ) that v(R,) = r mm . 

So far, we did not discuss what enters I \r) and thus H(r). In 
general, the letter F is used to describe a ratio between a radiative 
acceleration and gravity. The most common use is the so-called 
“electron gamma”. 


, _ «lhomO) — 0“ e L 

g(r) 4 k cm H G^ lon M» ’ 


( 22 ) 


which accounts only for the radiative pressure because of scat¬ 
tering of free electrons. Note that r 2 cancels out in the last ex¬ 
pression, since the acceleration because of Thomson scattering 
is defined as 


^thom(^*) • 


cr e L 


> ‘/.OllOf 


47TCWHr 2 
with the ionization parameter, 
n c (r) n e (r) 


q m n(r) = m H 


p(r) 


(23) 


(24) 


The mean atomic mass fR is constant in an atmosphere model. 
The two density factors can be replaced by the slowly varying 

mean particle mass p(r) — fR (l + pppj) such that 


flion i.b) — , 

P(r) 


1 

fR' 


(25) 


illustrating that q nm is nearly constant throughout the atmo¬ 
sphere. As a consequence, T e has only a very weak radial de¬ 
pendence^] 

While Eq. ( [22] ) is the common definition for F in the context 
of the Eddington limit, it is only a fraction of the total outward 
radiative acceleration a m d- Apart from the Thomson opacity cor¬ 
responding to scattering of free electrons, additional opacities 
originate in line (bound-bound) and continua (bound-free, free- 
free) transitions, contributing to the total radiative acceleration, 
a,ad is thus separated into three parts: 


^rad ^tlioni “b ^lines "b ^truecont- 


(26) 


The last term aWcont refers to the acceleration originat¬ 
ing from bound-free and free-free continuum transitions. As the 
electron opacity also forms a continuum, these two terms are 
combined and referred to as the “continuum”. In contrast, the 
term without electron scattering is sometimes called the “true 
continuum”. To avoid any confusion, we will adopt this notation 
in the following. 

Similar to Eq. ( |22| ), we can thus define a T, which corre¬ 
sponds to the total acceleration. 


r , flradM 1 4k 1 
r ra d (''):=—— = —-7T K v (r)H v (r) dv, 


g(r) g(r ) c p(r ) 


f K - 


(27) 


3 Instead of using < 7 ion , one also finds the specific electron-scattering 
coefficient ,v c = n e cr e /p in the literature (e.g., Mihalas 1978[l. 
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Table 1. Characteristics of different stellar atmosphere codes concerning their treatment of the (quasi-)hydrostatic regime 



PoWR 

CMFGEN 

FASTWIND 

TLUSTY 

radiative transfer 

comoving frame 

comoving frame 

CMF/Sobolev° 

static (obs. frame) 

blanketing 

full 

full 

approximative 

full 

temperature stratifica¬ 
tion obtained by 

radiative equilibrium* or 
thermal balance 

radiative equilibrium*’ 6 

thermal balance 

radiative 

equilibrium* 

photosphere 

quasi-hydrostatic 

quasi-hydrostatic 

quasi-hydrostatic 

hydrostatic 

(p, c)-update 

consistent 

injections* 

start iterations 6 

v = 0, p consistent 

radiative acceleration 
used in hydrostatic eq. 

fllll £Z ra d 

full Cly ad 

flcont approximation^ 

full a md 

Hurb parameter in hy¬ 
drostatic eq.? 

yes 

yes 

no 

yes 

domain connection 

continuous du/dr h 

0.75 tw/ 

0-1 iconic (adjustable) 

no wind domain 


Notes. <a) Comoving frame (CMF) method for main elements, Sobolev approximation for trace elements. Details are _ 

In PoWR and TLUSTY, the radiative equilibrium is considered in two different flavors, the so-called integral form, 4 nf k v (5 
and the flux consistency 4n f H v dv = <r SB T 


given in Puls et al. ( 2005 


(Hubeny & Lanz 


1995a 


■ J v ) dv = 0, 

Hamann & Grafener|2003). CMFGEN instead only uses the integral form 


l Hillier|2003 1 . <c) In CMFGEN, the thermal balance is ca lled “electron energy b alance” (EEB) and used to check convergence and superlevel 
assignments, but not explicitly for temperature corrections (Hillier & Miller 1998 Hillier 2003^. (< * After convergence of the statistical equations, 
the stratification in the quasi-hydrostatic part is adjusted and the model calculation is “restarted” from a gray temperature distribution. The total 
number of these “restarts”, which are also referred to as “injections”, must be specified beforehand. At the start of a model an iteration is 
done for the quasi-hydrostatic domain, where a Rosseland optical depth based on LTE opac ities is calculated along with the calculation of the 
temperature and density stratification. The details are explained in Santolaya-Rey et ah (19971. w The hydrostatic equation is part of the set of 
linearized equations, which are solved consistently in every iteration ( Hubeny STLanzj 1995b 1 A g) The continuum acceleration is approximated by 
a nonintegral term with a parameterized Rosseland opacity (Santolaya-Rey et a l.|1997p . 1)11 The continuous velocity gradient is the standard option 
in PoWR to find a connection point. Alternatively, the user can specify that the connection point is forced at v = f ■ u sonic . (l) Models before ~2013 
used v = 0.5 u son ic for the connection point (jMassey et al.|2013)l. 


where k v and H v are the opacity and Eddington flux at the fre¬ 
quency v, respectively, and c is the speed of light. This “full” 
T rad is then used in the quasi-hydrostatic domain. A similar ap¬ 
proach was also adopted by Lanz & Hubeny (20031 in their 
plane-parallel TLUSTY code, which is widely used for calcu¬ 
lating photospheric spectra of static atmospheres. 


2.3. The quasi-hydrostatic treatment in different stellar 
atmosphere codes 

A couple of non-LTE stellar atmopshere codes are used in the 
field of hot and massive stars. Widely used are CMFGEN ([Hillier 


|& Miller|1 998), FASTWIND ([Santolaya-Rey et aL| 1 997[>, 


in case of negligible stellar winds, TLUSTY (Hubeny & Lanz| 
|1995b| l. While these codes have a lot of similarities in their con¬ 
cepts, they also differ in certain approaches, such as their treat¬ 
ment of the radiative transfer or their superlevel approach. Re¬ 
cent comparisons between CMFGEN and FASTWIND (Massey 


In contrast to T e , which is already known at the start of a 
non-LTE model atmosphere calculation, with the exception of 
the exact value of <7j on (r), the full r iat j has to be calculated iter¬ 
atively. The radiative acceleration is calculated from the popu¬ 
lation numbers, which in turn depend on the radiation field and 
also on the density which itself is connected to F ra( ] via the hydro¬ 
static equation. To ensure a consistent solution, the velocity field 
in the quasi-hydrostatic domain is constantly updated as soon 
as the hydrostatic equation is violated by more than 5%. This 
turned out to be sufficient in all test cases, as long as T rad (r) < 1 
in the quasi-hydrostatic part. In reality, values larger than unity in 
the (quasi-)hydrostatic domain would reflect additional physical 
phenomena occurring in the photosphere (e.g., subphotospheric 
convection, Cantie llo et al,||2009) , which are not handled by 
our model atmospheres. The plane-parallel TLUSTY code treats 
only values of T ra( j < 0.9 in the hydrostatic equation ( [Lanz &| 
Huben y|2003| l as the hydrostatic equation fails for values larger 
than unity and numerical instabilities already occur when T ra< i 
gets close to this limit. We use a similar approach in PoWR for 
the quasi-hydrostatic equation, i.e., the value of T lat j is limited to 
0.9 for our velocity calculations. 


|et al.|2013] l for particular objects have shown a good agreement 
in the derived temperatures, but also revealed that the surface 
gravities derived with FASTWIND are abound 0.12dex lower 
than with CMFGEN. As the derived surface gravities are directly 
connected to the treatment of the quasi-hydrostatic layers, it is of 
major interest to understand the difference between the codes in 
this regime. 

Table [T] provides an overview of the features affecting the 
quasi-hydrostatic treatment of PoWR, CMFGEN, and FAST¬ 
WIND together with the completely hydrostatic treatment in 
TLUSTY. The latter provides an interesting test case in the limit 
of negligible mass-loss rates, which will be further discussed in 
Sect. [4] where we compare spectra from TLUSTY models with 
those from PoWR. In CMFGEN, TLUSTY model atmospheres 
can also be used as a start approach. 

The most striking difference in the quasi-hydrostatic treat¬ 
ment between the different wind codes (PoWR, CMFGEN, 
FASTWIND) is the update mechanism of the density and ve¬ 
locity stratification: 

- In PoWR, we update the velocity field as described in 

Sect.|2.2|as part of the main iteration as soon as the hydro- 
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static equation is violated by more than 5% or the specified 
value of r max is missed by more than a specified e T . For each 
update, F ra d(r) is applied in the hydrostatic equation, using 
flrad(f) from the comoving frame (CMF) calculations. With 
the new velocity field given, the stratification is updated ac¬ 
cording to the equation of continuity d3j. 

In CMFGEN ( jHillier & Miller|1998[ ), tne hydrostatic domain 
is first taken from a TLUSTY model or an old CMFGEN 
model, but can be updated with regards to the hydrostatic 
equation after a first convergence of the population numbers 
and the radiation field. After each such stratification update, 
sometimes referred to as “injection”, the model iteration has 
to be continued to reachieve convergence (see, e.g., |Martins| 
et al.|2012 for a brief description). Similar to PoWR, the ra¬ 


diative acceleration used in the hydrostatic equation is taken 
from the CMF ra diative transfer calculations. 

- In FASTWIND (Santolaya-Rey et al. |1997] >, the line con¬ 
tributions are instead neglected in the quasi-hydrostatic do¬ 
main and uses an approximate handling for the continuum 
acceleration, including both Thomson and true continuum 
term. This simplification allows FASTWIND to set up the 
quasi-hydrostatic domain in an iteration right at the start of a 
model, in line with the goal of the code of being significantly 
faster than PoWR and CMFGEN. 


Furthermore, all three wind codes use a different criterion for 
the connection point between the quasi-hydrostatic and the wind 
domain. While FASTWIND and CMFGEN use a fixed fraction 
of the sound speed to define this connection point, PoWR puts 
the condition of a smooth velocity gradient between the two 
regimes. 

Instead of treating only the quasi-hydrostatic part consis¬ 
tent, one could also aim at a hydrodynamic ally consistent treat¬ 
ment throughout the whole model atmosphere, i.e., including 
the wind domain. These calculations are typically not included 
in state-of-the-art stellar atmosphere codes, especially as a pre¬ 
scribed /Maw provides a good approximation for the outer wind 
regime of OB stars. A self-consistent calculation of the ve¬ 
locity field in the whole stellar atmosphere has been imple¬ 
mented as a nonstandard option in PoWR in order to discuss 
WR winds ( jGrafener & Hamann||2005| |2008| ). However, such 
models require much longer computation times and this option 
is not required for the purpose, in which we focus on the quasi¬ 
hydrostatic regime. The implementation presented here is com¬ 
pletely sufficient in achieving a self-consistent stratification of 
the quasi-hydrostatic layers. Nevertheless, a revised implemen¬ 
tation of completely hydrodynamically consistent model atmo¬ 
spheres in the PoWR code will be addressed in detail in a future 
work. 


2.4. The effective gravity 

The hydrostatic Eq. ([5]) can be written in an even shorter form 
with the definition of the effective gravity, 

GM , 

geff(/) := — 2~ [1 - F(r)]. (28) 

r i 

To avoid any confusion, we will use g grav (r) instead of g(r) with 
g* = ggrav C^*) to denote the full gravity from here on. The effec¬ 
tive gravity g c n is a fundamental fitting parameter for O- and B- 
stars when comparing models with observations. The pure grav¬ 
itational acceleration g gmv (r) = GM,r 2 , from which the stellar 
mass is derived, can only be calculated if r rac j is known. There¬ 
fore, g grav is model dependent. In these models an accurate treat¬ 
ment of the complete radiative pressure is required. Correcting 
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the observed g e ff for Thomson pressure r e only neglects a large 
part of the radiative pressure, and leads to an underestimated stel¬ 
lar mass, as we will demonstrate in Sect. [3] In the cases of neg¬ 
ligible radiative pressures only do log g grav and log g e ff become 
indistinguishable. Otherwise, spectroscopic masses are underes¬ 
timated if calculated directly via K; g e ff/G. 

As the shape of a spectral line like He) depends on g e ff and 
not g grav , there is a parameter degeneracy. A star with a weaker 
radiative pressure will produce the same line profiles as a star 
with a higher radiative pressure and a higher gravity. 

One of the major reasons for a strong model dependency of 
the stellar mass is the large impact of the opacity on the radiative 
pressure. It is significant which opacities are taken into account 
in the model. Neglecting elements such as Fe or Ne will result 
in a lower r ial j and is not legitimate, even if the corresponding 
element is not under consideration in the intended analyses. The 
inferred stellar mass therefore depends not only on the adopted 
stellar parameters, but also on the details of the model atoms. 

The adopted Doppler broadening velocity i>dop used in the 
comoving-frame calculations also has a significant impact: 
larger velocities allow atoms to absorb photons in a wider fre¬ 
quency range, and thus generally increase the radiative pressure 
in the atmosphere. 

Even though g s ff(r) and r rad (r) are calculated for each depth 
point in the PoWR code to fulfill the hydrostatic equation, there 
is a need for a certain “reference value” , which can be used 
for comparisons. As the emergent spectrum is mainly formed at 
tross ~ 2/3 in the photosphere, the first idea would be to use 
geft at this value. However, such a fixed point might already be 
located in the wind. Therefore we define a weighted mean of F ra( j 
as 

T sonic 

r ra d := J' r rad (r Ross ) e~ TRoss dr Ross . (29) 

^"max 

The upper limit of the integral r son ; c usually denotes the optical 
depth of the sonic point. However, we do not allow that this value 
drops below 0.1, even in the rare cases where the actual sonic 
point would be in such a regime. The calculated value of I ' rad is 
used to relate logg eff and log g, = logg grav (R„) via 

log geff = log g* + log (l -F rad ). (30) 

The radial dependence of r rad for a supergiant test model is 
illustrated in Fig. [I] where it is plotted against the total Rosse- 
land opacity. It is evident that in the quasi-hydrostatic regime 
not only the line transitions, but also the atomic continuum tran¬ 
sitions cannot be neglected, while this quickly changes as soon 
as we enter the supersonic domain, where F cont consists only of 
r e . One can see that the calculated value of I ' rad is indeed rep¬ 
resentative for the difference between g grav and g e ff in the photo¬ 
sphere. Values of F rad for a series of test models are shown and 
discussed in Sect. [3] 

The reference radius for all values of g e ff given in this work 
is R *. However, other radii are also used in the literature, typi¬ 
cally Ri — R (r Ross = |) or the radius of the sonic point R s . The 
latter roughly indicates the outer end of the quasi-hydrostatic do¬ 
main. It is therefore interesting to check how much the values of 
geff and g grav are affected when referring to different radii. As 
this effect is largest for supergiants, we examine an O-star model 
with T* = 32.5 kK, log L/L Q = 5.6, logg grav (R*) = 3.25 and 
logM = -5.75 [M 0 /yr]. For this model, we have Ri = 1.06R* 
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Fig. 1. Radial dependence of T ra d (red, solid line, upper panel) for our 
supergiant test model, plotted over the Rosseland optical depth r Ross . 
The mean value obtained with Eq. {29} is indicated as a dashed-dotted 
line. For comparison, T e and the total continuum contribution T cont are 
also shown. In the lower panel, the difference between g grav and g c ff is 
illustrated: the red solid curve marks the actual difference at each optical 
depth, while the dashed-dotted line denotes the difference using T rad . 
The gray dashed vertical line denotes the location of the sonic point. 


and the sonic point is located at R s = 1.1 I R„, leading to dif¬ 
ferences of 0.05 dex and 0.09 dex in log g. Especially the latter 
value is approximately what can be achieved by accurate fitting 
of high quality spectra and, therefore, we have to take the ref¬ 
erence radius into account when analyzing these objects. This 
effect is lower for giants and dwarfs, namely on the order of 
0.03 dex and 0.01 dex for the sonic point, respectively. 

In the literature, the term g e ff may not always have the same 
meaning. Especially when dealing with observations, it is often 
not clearly stated whether values termed as log g have been cor¬ 
rected for radiative acceleration. Values that have been corrected 
for centrifugal acceleration in a statistical sense are sometimes 
labeled as “true” gravities g lruc or g c with the relation 


gc =£ + 


(v sin i) 2 

R , 


(31) 


3. Results and discussion 

3. 1. Test model details 


To illustrate the impact of accounting for the full radiative pres¬ 
sure consistently, we calculate a set of PoWR models for a fixed 
temperature of T, = 32.5 kK, which corresponds to a late O-star. 
For this temperature, we calculate three types of models with dif¬ 
ferent surface gravities, corresponding to the luminosity classes 
I (supergiant), III (giant), and V (dwarf). As we focus on the 
quasi-hydrostatic part, we adopt only schematic parameters for 
the wind, i.e., all models have the same terminal wind velocity 
of v m = 2000 km/s. The outer model boundary is set to 100 R t . 
The velocity field in the wind domain can be prescribed by a 
f5 —law with fi = 0.8 (see Eq.[2}, which is accurate enough for our 
purposes as we do not want to analyze the outer wind. The lu¬ 
minosities are “typical” representatives of their class according 
to Martins et al. ( |2005} . We chose the surface gravities simi¬ 
larly, but with slight adjustments of up to 0.1 dex, correspond¬ 
ing to the closest grid point in the TLUSTY O-star model grid 
( |Lanz & Hubeny|2003[ l for a later comparison (see Sect. [4}. We 
adopted mass-loss rates from Vink et aT] {2000}, and the wind 


is assumed to be smooth (no density contrast, i.e., D = 1 in the 
notation of Hamann & Koesterk e|l998} . The chemical composi¬ 
tions are taken to be solar. Here, we use the abundances inferred 
Grevesse & Sauval ( 1998} instead of the newer ones obtained 


by 
by 

the TLUSTY O-grid. All parameters of the models are compiled 
in Table|2] 


Asplund et al. (2009) to have identical values to those used in 



io gl0 y?/R„ -1) 


Fig. 2. Acceleration stratification for a model with a consistent quasi¬ 
hydrostatic domain (Type (b) models). The wind acceleration (thick red 
diamond line) is compared to the repulsive sum of inertia and gravita¬ 
tional acceleration g(r) (black line). All terms have been normalized to 
g(r). 


( |Repolusr et al. 2004). The term labeled as g in this equation 
is sometimes called g e ft (e.g., |Massey et al,||2013} . We would 
like to stress that this is not identical to our definition of the 
effective gravity in this work. The variable g in {31} refers to 
the gravitational acceleration specified in a non-rotating stellar 
atmosphere model, which has not been reduced by radiative ac¬ 
celeration. As such, g would be referred to as g gl - av , following the 
notation of this work. Caution is therefore adviced when dealing 
with the term “effective gravity”, as its definition might differ 
significantly between different authors. 


For each of the three luminosity classes, we calculate three 
models, which we refer to as (a), (b), and (c) in the following, 
making nine models in total period. Models (a) and (b) have 
a fixed g gmv , which corresponds to the respective luminosity 
class and differ only in the treatment of the radiative accelera¬ 
tion in the quasi-hydrostatic part: models (a) are only calculated 
with r e , while models (b) include the full r ra j, i.e., they ac¬ 
count for the complete radiative pressure. Typical examples are 
shown in Fig. [2] and Fig. [3] where we plot the acceleration strat¬ 
ifications for a dwarf (class V) model. One can see that Fig. [2] 
shows the consistent (b)-model, where the sum of the total ra- 
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MA] 


Fig. 4. Balmer line profiles of H(5 (left) to Ho- (right) for test models with fixed g grav . The spectra from the (a)-type models with a quasi-hydrostatic 
stratification calculated using only the Thomson term T e are shown as blue dashed line. The red solid lines display the profiles from the (b)-type 
models accounting for the full radiative term r ra d. 



log l0 (R/R, - 1) 


Fig. 3. Same as Fig.[2] except with a velocity law in the quasi¬ 
hydrostatic domain accounting only for T c (Type (a) models). While 
the wind domain is not affected, the acceleration balance in the quasi¬ 
hydrostatic part is significantly different. 


diative acceleration and gas pressure equals gravitation in the 
quasi-hydrostatic domain. The result of the (a)-model is shown 
in Fig. [3] Now only I e and gas pressure are taken into account for 
the hydrostatic equation and thus the sum of the total radiative 
acceleration and gas pressure is larger than the local gravity in 
the subsonic part. The outer parts of both models are very similar 
as both use the same prescribed /3-velocity law. 

Article number, page 8 of 


In models (c), like in models (a), we account only for the 
Thomson term T e in the quasi-hydrostatic domain, but instead 
of specifying logg grav , we specify the effective gravity logged 
and adopt its value from the corresponding model (b), calculated 
from Eq. ( |30| . An overview of the basic similarities and differ¬ 
ences between all three model types is given in Table[3] Since 
models (b) and (c) of a given luminosity class have the same 
effective gravity, the wings of pressure-broadened lines are ex¬ 
pected to be identical in both models. However, because of dif¬ 
ferent treatment of the quasi-hydrostatic domain, the actual sur¬ 
face gravities g gmv and spectroscopic masses M, implied from 
both models will differ. 


3.2. Comparison with fixed g grav 

Figure[4] illustrates the impact of accounting for the full radia¬ 
tive pressure on prominent Balmer lines (left to right: Hd, Hy, 
H/3, Her) for dwarfs (upper panels), giants (middle panels), and 
supergiants (lower panels). The impact of either including only 
the Thomson term F e (model a, blue dashed lines) or includ¬ 
ing the full radiative term F ra( j (model b, red solid lines) can be 
seen in the line wings. Because of their larger outward pressures, 
the quasi-hydrostatic domains of models (b) are less dense than 
those of models (a), and, as a consequence, the line wings ob¬ 
tained by models (b) are narrower. For the same reason, the value 
of geff is always smaller in models (b) compared to models (a) 
(see below). 

A further inspection of Fig.[4]reveals that the difference be¬ 
tween models (a) and (b) increases with luminosity class. A look 
at Eq. (27 1 reveals that lower values of log g grav will increase 
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Table 2. O-star test model parameters 


Fuminosity class 

i 

III 

V 

T, [kK] 


32.5 


loggeff [eras 4 ]” 

2.47 

3.02 

3.79 

logggrav [ems 4 ]" 

3.25 

3.50 

4.00 

T ra d 

0.83 

0.67 

0.39 

T e 

0.40 

0.22 

0.07 

R * [*o] 

20.0 

14.1 

8.1 

log M [Mq yr" 1 ] 

-5.75 

-6.25 

-7.1 

log L[Lq] 

5.60 

5.30 

4.82 

Mt [Mq] 

25.9 

23.0 

24.1 

p 


0.8 


Uoo [km/s] 


2000 


Uturb [km/s] 


10 


u dop [km/s] 


20 




0.704 




0.282 


V 

2.78 x 10" 

3 

V 

8.14 x 10" 

4 

V 

7.56 x 10" 

3 

XMg b 

6.45 x 10" 

4 

x M b 

5.56 x 10" 

5 

x Si b 

6.96 x 10“ 

4 

V 

6.10 x 10"' 

6 

Xs b 

4.79 x 10“ 

4 

V b,c 
^Fe 

1.34 x 10" 

3 


Notes. {a> Values refer to the consistent models. Comparison model re¬ 
sults are given in Tablejd ] (b) Solar abundances, as obtained by Grevesse 


& Sauval 

kO 

kO 

OO 

i, specified here as mass fractions. (t) Fe includes also 

the further iron 

group elements Sc, Ti, V, Cr, Mn, Co, and Ni. See 

Grafener et al. (2002) for relative abundances. 


Table 3. O-star test model type overview 


Model type (a) (b) (c) 


8 grav 

fixed 2 

fixed 2 

calculated 

r ; 

T e 

T ra d 

r e 

geff 

calculated 

calculated 2 

fixed 2 


Notes. (,) T value used to relate g geff and g grav <2) Type (a) and (b) models 
within the same luminosity class have the same g grav (J| A type (c) model 
has the same g geff as the type (b) model of the same luminosity class 

r ra d, which is indeed much larger for the supergiant than for the 
dwarf. Even though the difference between r ra d and F c is not 
getting much stronger with lower log g„ rav , as T e also changes 
proportional to both values are significantly larger for the 
supergiant and thus the difference in the derived g cfr value in¬ 
creases, which is reflected in the increasing line wing difference 
in Fig. [4] 

Even for our dwarf models (log g grav = 4.0, upper panels), 
where the dashed and solid lines can only barely be distin¬ 
guished, we obtain a difference of A logged = 0.21 dex when 
accounting for the full radiative pressure via r ra d (logged = 3.97 
for model (a) vs. logged = 3.78 for model (b)). For our gi¬ 
ant models (logg grav = 3.5, middle panels) the difference in¬ 
creases to 0.37 dex (3.39 vs. 3.02). Finally, for the supergiant 
models (logg glav = 3.25, lower panels), a formidable difference 
of 0.56 dex (3.03 vs. 2.47) is obtained. Even though the impact 
on the spectral appearance of the hydrogen lines might not be 


particularly striking in Fig. [4] the differences in the correspond¬ 
ing values of log g e g are quite remarkable. 



x/A x/A 


Fig. 5. The Hy and H/f lines (left and right panels, respectively), as 
obtained from models (b) and (c) (red solid and green dotted lines, re¬ 
spectively), both with an effective gravity of logg eff = 2.47. Although 
the strengh of the line wings is the same, the two models imply masses 
and gravities, which strongly differ because of the different treatment of 
the radiative pressure (cf. Table|4j 


3.3. Comparison with fixed g e tf 

In the models (a) and (b) discussed above, we specify g grav , 
T„, and L, and therefore the radii, the spectroscopic masses are 
known a priori. When dealing with observations, however, it is 
getf, and not g grav , which is empirically measurable. Thus, one 
would fix the effective gravity that best reproduces the observa¬ 
tion, and let the actual gravity g grav be a model output. Based 
on the detailed physics implemented in the model atmosphere, 
the gravity g grav and the spectroscopic mass M, follow from 
loggeff- This scenario is reflected in models (b) and (c), which 
are calculated with identical effective gravities, but with a differ¬ 
ent treatment of the radiation pressure in the quasi-hydrostatic 
domain. Fig. [5] shows the Balmer lines Hy (left panel) and H/j 
(right panel) as obtained from the supergiant models (b)„ which 
account for r ra d (red solid line), and (c), which only include F e 
(green dotted line). As both models have the same g e ff, the line 
wings obtained from both models can hardly be distinguished. 
However, some differences are seen in the line cores and in the 
helium and metal lines. 

To illustrate the effect of this different treatment, we show the 
velocity stratification for all three supergiant models in Fig. 6] 
All models agree in the outer part where we have prescribed the 
/1-law. However, in the inner part it becomes evident that fixing 
geff indeed leads to relatively similar velocity stratifications with 
both the consistent r ra d-approach and the r e -approach, which is 
reflected in the agreement of line profiles between the models 
(b) and (c). If one fixes instead the more fundamental parameter 
g grav , the resulting velocity stratification is significantly differ¬ 
ent in the inner part, if accounting only for T e , including a huge 
difference in the radial location of the sonic point. These dis¬ 
crepancies explain the spectral differences between models (a) 
and (b) shown in Fig. [4] 

Despite the similarities between the models (b) and (c) seen 
in the line wings and in the velocity law, their values of log g grav 
and M, differ significantly. Tableffl shows the values of logg grav 
and M* for the dwarf, giant, andsupergiant models. As could 
be anticipated from Eq. ( [30] ). the models that account for the full 
radiative pressure (models b) have larger values for log g grav and 
M t because their T, which is used to relate g e ff and g grav , is larger 
in these models. The spectroscopic masses deviate by roughly a 
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log (r/R, - 1) 

Fig. 6. Velocity stratification for three supergiant test models: the con¬ 
sistent quasi-hydrostatic model (b) (red solid) is compared to a model 
(a) (blue dashed) with the same g grav , but where only T e is taken into 
account in the quasi-hydrostatic part (blue). The third model (c) (green 
dotted) shares g eff with model (b), except it accounts only for T e and 
thus has a different g grav . 

Table 4. Deduced stellar masses with different quasi-hydrostatic ap¬ 
proaches 


Luminosity class 

I 

III 

V 

logoff [cms' 2 ] fl 

2.47 

3.02 

3.79 

logg g mv(r e )[cms" 2 ] 6 

3.00 

3.25 

3.82 

logg grav (r rad ) [cm s' 2 ] 

3.25 

3.50 

4.00 

M.(Te) [M Q ] b 

14.5 

12.9 

16.2 

4Tt(T rad ) [Mg] 

25.9 

23.0 

24.1 


2.4 1 is fixed for all models here 


Notes. (a) The effective gravity (cf. Sect. 

(il These values are derived from models that were calculated as if only 
the Thomson radiative pressure would enter the hydrostatic equation. 


marginally, we can deduce from Eq. ( [22] ) that the product of T e 
and M, will be approximately the same for both models. Hence 
we get an expression that allows us to replace T^f with the known 
i.e., 


w 


r w sb M * r w 


(c) 


M{ 

Using that in ([32| yields 


Mi e) - M[ b) r® » Mi b) 


(l - Trad) 


Mi c) « Mi h) 


( r rad - T?') 


(33) 


(34) 

(35) 


This means that masses of the (c) models only depend on the 
difference between T la d and T e in the consistent (b) models, not 
on the absolute values. The difference is comparable for all three 
luminosity classes, and so is the mass deviation in Table[4] 



li A li A 


Fig. 7. Left panel: the spectrum around Hy for the hydrostatically con¬ 
sistent supergiant test model (red solid line) is compared to a simi¬ 
lar PoWR model calculated without iron group elements (black dot¬ 
ted line). Right panel: the same test model is calculated with a higher 
Doppler broadening velocity Udop (green dashed line) and compared to 
the original. 


factor of two. Interestingly, there is no clear trend for increasing 
deviation with increasing luminosity class. Intuitively, one might 
expect that the deviation would have to be larger for the super¬ 
giant model, compared to the dwarf model, as we obtained it in 
spectral comparison of the models (a) and (b) which were shown 
in Fig. [4] However, the calculation of the (c) models is quantita¬ 
tively different from the (a) models. Even though both consider 
only T e for obtaining the velocity field in the quasi-hydrostatic 
part, the (c) models are specifically designed to reproduce the 
geff-value obtained with the full T ra( j , while the (a) models have 
a completely different approach, with the ggrav-value being iden¬ 
tical to the (b) models. In fact, we can estimate the masses of the 
(c) models with an easy calculation. Starting from the require¬ 
ment that both, (b) and (c) models should have the same g e ff, it 
immediately follows via Eq. ( |28j ) that 

M< c) (l-rM) = M®(l-f®). (32) 


The small superscripts indicate the value from the corresponding 
model family, i.e., M; is short for M*(T e ) in the (c) models. As 
we do not know Te in advance, since it contains M*(T e ) itself 
by definition (22 1 , we need to replace it with a value from the 
(b) model. Because the luminosity L is the same in the (b) and 
the (c) models and the ionization parameter < 7 j on only changes 
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3.4. Blanketing and Doppler velocity influence 

It is unfortunate that the spectroscopic mass greatly depends on 
various parameters adopted in the calculation. Two important ex¬ 
amples will be discussed here. The first is the importance of the 
iron group elements. Since they are a dominant source for opac¬ 
ity in the atmosphere of a massive star, their abundances signifi¬ 
cantly affect the radiative pressure, and thus the inferred spectro¬ 
scopic mass. Therefore iron group elements have to be included 
in model atmosphere calculations, even if these models are only 
used to analyse spectral regimes without visible iron lines. The 
effect, which is called “line blanketing”, does not only affect the 
temperature stratification, but also the density structure of a stel¬ 
lar atmosphere. This is illustrated in left panel of Fig. [7] where 
we compare the Hy line from our supergiant model (see Table[2] 
for parameters) to that of an identical model without iron. It im¬ 
pressively demonstrates that neglecting elements, which greatly 
contribute to the total opacity, is not legitimate. When the iron 
elements are not included, the radiative pressure is smaller, and 
the line wings appear broader in the synthetic spectrum, lead¬ 
ing to an underestimation of the gravity by ~ 0.3 dex. These ele¬ 
ments should thus be included in any consistent calculation, even 
if there is no observational abundance indicator. Using “typical” 
abundances from the local region or similar objects will usually 
cause a smaller error than neglecting such elements completely. 
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Another quantity that is of major influence on the spectral ap¬ 
pearance of OB star atmosphere models is the adopted Doppler 
broadening velocity i^op- This velocity is used in the comov¬ 
ing frame calculations and reflects the combined influence of 
the thermal and micro turbulent velocities. While the thermal ve¬ 
locity is calculated for each element as a depth-dependent man¬ 
ner in the formal integral, the comoving frame calculations are 
currently using a constant Udop for all elements. For the spectral 
appearance of Wolf-Rayet atmospheres the value of i;j op is of 
minor importance, but it has a notable impact on the spectra of 
our OB models, as it affects the total radiative pressure in the 
atmosphere. This is illustrated in the right panel of Fig. [7] where 
we again compare the region around Hy for our supergiant test 
model, but now compare with a model that has been calculated 
with a larger Doppler velocity. In the O-star regime, the value 
°f fdop = 20km/s has proved to be sufficient. This demonstrates 
that it is imperative to choose a i>dop reflecting the true thermal 
and turbulent velocity in the stellar atmosphere, despite the sig¬ 
nificantly longer computing times in the comoving frame calcu¬ 
lations that come along with smaller values of i^op- 


sence of emission lines is not sufficient to deduce that mass loss 
is negligible. 



nk 



x/k 


Fig. 8. Profiles of the Her and H/3 lines from PoWR models with differ¬ 
ent M. In the limit of low mass-loss rates, the wind effect on the spec¬ 
trum becomes negligible. For an easier comparison of the line wings, 
the electron redistribution (see Fig.|9]f is switched off in these simula¬ 
tions. 


4. Comparison with TLUSTY 


The TLUSTY code ( |Hubeny & Lanz||1995b| l provides plane- 
parallel model atmospheres. The synthetic spectra that can be 
calculated from these models, e.g., with the SYNSPEC program 
from the same authors, are widely used for the analysis of pho- 
tospheric spectra of OB-type stars. In the following section, we 
compare PoWR and TLUSTY model results. In all cases where 
spectra based on TLUSTY models are shown, they were ob¬ 
tained with the SYNSPEC code and are labeled as TLUSTY. 
0 


While both PoWR and TLUSTY are non-LTE codes, a ma¬ 
jor difference is that the PoWR models include the stellar wind. 
PoWR therefore adopts a spherical geometry, while TLUSTY 
assumes a plane-parallel geometry. As thoroughly discussed by 
|Lanz & Hubeny| ( 2003[ ), the assumption of static, plane-parallel 
atmospheres is a fairly solid approximation for the photospheres 
of OB stars, which show negligible signs of stellar winds and 
curvature effects. One can therefore expect that in the limit of 
negligible mass-loss rates and small scale heights (see Eq.|4|, the 
PoWR spectra would be close to the TLUSTY results. To check 
that PoWR models of OB-type stars with negligible winds show 
a good agreement with the corresponding TLUSTY models, we 
dedicate this section to a comparison of dwarf, giant, and su¬ 
pergiant PoWR models with their TLUSTY counterparts, focus¬ 
ing on the pressure-broadened Balmer lines. In all cases, PoWR 
models with the consistently treated quasi-hydrostatic domain 
are used, i.e., those models that we referred to as type (b) above. 


Since we calcuated the models used in Sect.[3]with “typi¬ 
cal” O-star mass-loss rates, they are not expected to provide a 
very close agreement with the TLUSTY models. The discrep¬ 
ancies are most prominent in the Her line for the supergiant and 
giant models, while small differences are also found in the other 
Balmer lines. In Fig. [8] we illustrate the effect of reducing the 
mass-loss rate. It is evident that the spectral appearance con¬ 
verges in the limit of small mass loss, and it is at this limit where 
we expect to obtain the closest agreement with TLUSTY mod¬ 
els. Furthermore, the sequence of models illustrates that the ab- 


4 The models from the TLUSTY O- and B-star grids (including their 
SYNSPEC spectra) can be obtained from the TLUSTY website at 
http://nova.astro.umd.edu/ 
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Fig. 9. Line profiles of Hy and Hb for the same supergiant PoWR model. 
The red line denotes the spectrum that is obtained when switching off 
the electron redistribution in the formal integral, while the black line 
indicates the “normal” output, which includes the redistribution of free 
electrons. 


Another PoWR feature, which would hinder a comparison 
with the TLUSTY models, is frequency redistribution by Thom¬ 
son scattering ( jMihalas||1978| >. Given their large thermal veloc¬ 
ities (-1000 km s 1 ), free electrons can scatter photons to sig¬ 
nificantly different wavelengths. Hence, photons that are trapped 
in an optically thick line core can be Doppler-shifted to the line 
wing or adjacent continuum, from which they can freely emerge 
and become visible as an excess emission. Especially at lower 
surface gravities and in spectral domains with a high density of 
spectral lines, the effect of frequency redistribution often results 
in a noticeable pseudocontinuum. When the synthetic flux is nor¬ 
malized relative to the continuum, this has the appearance of a 
continuum offset in the normalized spectrum. To illustrate this 
effect. Fig. [9] shows the Hb and Hy lines (left and right panels, 
respectively) for the consistent supergiant model including fre¬ 
quency redistribution (black solid line) and without (red solid 
line). Since TLUSTY does not account for this effect, we dis¬ 
able it for the sake of comparison. However, we stress that the 
importance of this effect has been demonstrated in various stud- 


ie; 


(Hummer & Mihalas 1967 Auer & Mihalas 1968). 


5 It may seem arbitrary to rectify the synthetic flux using the contin¬ 
uum before accounting for the redistribution. However, the redistributed 
flux varies significantly around spectral lines, while the “unredistributed 
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Fig. 10. Balmer line comparison between the PoWR models (red. solid line) and the TLUSTY models (black, dashed line) for dwarfs (upper 
panels), giants (middle panels), and supergiants (lower panels). The stellar parameters used in the PoWR models are as in Table[2] but with 
negligible mass-loss rates (log M = -9 [M 0 yr' 1 ]), small terminal velocities (500km s _1 ), and large radii, to avoid wind and curvature effects. 


The upper, middle, and bottom panels of Fig.[l0]show a com¬ 
parison between the TLUSTY (black, dashed line) and PoWR 
(red, solid line) for the dwarf, giant, and supergiant models, re¬ 
spectively, where we compare the first four Balmer members 
(from left to right: Hd, Hy, H/j, Ha). The PoWR model parame¬ 
ters are identical to those compiled in Tablepl but with a negli¬ 
gible mass-loss rate of log M = -9.0 [M Q yr^ 1 ] and a small ter¬ 
minal velocity of 500 km s '. Furthermore, the radii of all stars 
were set to large values to diminish possible curvature effects 
and to approach a plane-parallel geometry. 

It is evident that a very good agreement in the line wings is 
obtained between the PoWR and TLUSTY models for all three 
luminosity classes. In particular, the wings of the Balmer lines 
are barely distinguishable. Hence, the application of either mod¬ 
els would result in practically the same surface gravity, provided 
that the mass loss is really negligible. Interestingly, pressure 
broadening in the TLUSTY O-grid models from |Lanz & Hubeny 
( |2003|> uses analytical a pproximations to numerical calculations 
(cf. Hubeny et al.|1994] appendix B), while PoWR performs in¬ 
terpolation over pressure broadening tables ( |Lemke 1997 1 for 
hydrogen lines. Given the good agreement, both methods seem 
to be adequate for the studied parameter regime, i.e., nondegen¬ 
erated stars. 

While a comprehensive comparison between PoWR and 
TLUSTY along all spectral lines is beyond the scope of the cur¬ 
rent paper, we note that there are significant differences regard- 


continuum'’ is a slowly varying function of A and therefore much more 
appropriate for normalization. Moreover, electron scattering contains 
vital information regarding the physics in the stellar atmosphere, and 
should not be removed by normalization. 
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ing helium and metal lines. Not only does their strength differ by 
up to a factor of two between TLUSTY and PoWR, but some¬ 
times - as visible in the supergiant comparison - they may even 
appear in emission from one code, and absorption from the other. 
Because of the non-LTE conditions, even small differences in the 
stratification can have large effects on such weak lines, e.g., be¬ 
tween two PoWR models differing only slightly in their tempera¬ 
ture or gravity. Indeed, we find similar discrepancies in the small 
lines. This can be seen in Fig.[i] where we compare the F ra( j- 
and F e -model spectra. Therefore one has to be careful deducing 
parameters, e.g., abundances, from only one particular line. In¬ 
stead, several lines of more than one ionization stage should be 
compared if available. 


5. Summary and conclusions 

We present a set of stellar atmosphere models calculated with 
the most recent version of the PoWR code, using different ap¬ 
proaches for the quasi-hydrostatic regime. In the limit of small 
mass-loss rates, we also compared the PoWR models to the 
plane-parallel TLUSTY atmospheres. 

We conclude that a proper treatment of the quasi-hydrostatic 
regime is imperative for OB-type star modeling. For a consis¬ 
tent solution, the full radiative acceleration has to be taken into 
account, including line and continuum contributions. 

The spectroscopic masses will be severely underestimated, 
by a factor of roughly 2, if models are used that account only for 
the radiation pressure on free electrons in the quasi-hydrostatic 
domain. This holds for all luminosity classes. 
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Omitting elements that significantly contribute to the to¬ 
tal opacity compromises the density stratification in the quasi¬ 
hydrostatic domain, and consequently leads to an inconsistent 
stellar masses. Specifically neglecting important elements, such 
as Fe, in the models is by no means legitimate. 

The effect of mass loss on the spectra of typical OB stars may 
seem subtle, but even small mass-loss rates can change the line 
profile and thus might affect deduced gravities. The absence of 
emission lines in an observation does not imply negligible mass 
loss. 

In the limit of small mass-loss rates and vanishing curvature 
effects, the emergent spectra of the PoWR model atmospheres 
generally agree very well with TLUSTY models calculated with 
the same stellar parameters. 
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